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ABSTRACT 

We explore the formation of dusty gas walls induced by a circumnuclear starburst around an active 
galactic nucleus (AGN). We concentrate our attention on the role of the radiation force by a starburst 
as well as an AGN, where the effects of optical depth of dusty gas are taken into consideration. First, we 
solve the hydrostatic equations in spherical symmetry coupled with the frequency-dependent radiative 
processes, to demonstrate that a geometrically thin, optically thick wall forms due to the radiation 
pressure by a circumnuclear starburst. Next, in two-dimensional axisymmetric space, we analyze the 
configuration and the stability of geometrically thin walls which are in balance between radiation pressure 
and gravity. As a result, it is shown that the radiation force by the circumnuclear starburst works to 
stabilize optically thick walls surrounding the nucleus. In the case of a brighter starburst with a fainter 
AGN (e.g. L SB /M SB >10[L Q /M Q ] and Lagn <lO n L ), there form double walls, an inner one of which 
is located between the nucleus and the circumnuclear starburst, and an outer one of which enshrouds 
both the starburst regions and the nucleus. The total extinction of both walls turns out to be larger 
for a brighter starburst, which is Ay ~ 10 for -Lsb/A^sb >10 2 [£q/M q ]. As a consequence, double walls 
could heavily obscure the nucleus to make the AGN type 2. The outer wall may provide an explanation 
for the recent indications for large-scale obscuring materials in Seyfert 2's. Also, it is predicted that the 
AGN type is time-dependent according to the stellar evolution in the starburst, which shifts from type 
2 to type 1 in several times 10 7 yr owing to the disappearance of walls. In contrast, if the AGN itself 
is much brighter than the starburst as a quasar is, then neither wall forms regardless of the starburst 
activity and the nucleus is likely to be identified as type 1. To conclude, the radiatively-supported gas 
walls could be responsible for the putative correlation between AGN type and the starbursts, whereby 
Seyfert 2 galaxies are more frequently associated with circumnuclear starbursts than type 1, whereas 
quasars are mostly observed as type 1 regardless of star-forming activity in the host galaxies. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: nuclei — galaxies: starburst — 
quasars: general — radiative transfer 



1. INTRODUCTION 



. There has been a good deal of evidence on an obscuring 
$— i torus of subparsec scale, which surrounds an active galac- 
fic nuclei (AGN) (Antonucci 1984; Wilson, Ward, & Haniff 
" 1988; Barthel 1989; Blanco, Ward, & Wright 1990; Miller 
& Goodrich 1990; Awaki et al. 1991; Storchi-Bergmann, 
Mulchaey, & Wilson 1992). The obscuring torus is thought 
to be responsible for the dichotomy of AGN type in the 
context of the unified model (Antonucci 1993, for a re- 
view). However, the origin and physical structure of the 
torus has not been well elucidated, although some intrigu- 
ing models are proposed by Krolik & Begelman (1988), 
Pier & Krolik (1992a, 1992b, 1993), Efstathiou & Rowan- 
Robinson (1995), and Manske, Henning, & Men'shchikov 
(1998). 

Recent observations on AGN hosts have gradually re- 
vealed that the properties of host galaxies of Seyferts are 
intrinsically unlike between type 1 and type 2 (Heckman 
et al. 1989; Maiolino et al. 1995, 1997, 1998a; Perez- 
Olea & Colina 1996; Hunt et al. 1997; Malkan, Gorjian, 
& Tarn 1998; Storchi-Bergmann, Schmitt, & Fernandes 
1999; Storchi-Bergmann et al. 2000; Gonzalez Delgado, 
Heckman, & Leitherer 2001). The Seyfert 2 galaxies are 
more frequently associated with the circumnuclear star- 



bursts than type l's. In contrast, quasars (QSOs) are 
mostly observed as type 1, although the QSO hosts often 
exhibit vigorous star-formation activity (Barvainis, An- 
tonucci, & Coleman 1992; Ohta et al. 1996; Omont et 
al. 1996; Schinnerer, Eckart, & Tacconi 1998; Brotherton 
et al. 1999; Canalizo & Stockton 2000a, 2000b; Dietrich 
& Wilhelm-Erkens 2000). [It is noted that the QSO hosts 
are fainter than QSO nuclei themselves (McLeod & Rieke 
1995b; Bahcall et al. 1997; Hooper, Impey, & Foltz 1997; 
Crawford et al. 1999; Kirhakos et al. 1999; McLure et al. 
1999; McLure, Dunlop, & Kukula 2000).] Such correlation 
with circumnuclear starbursts seems beyond understand- 
ing based upon the picture of the unified model where the 
bifurcation of AGN type is simply accounted for the ori- 
entation of the nucleus surrounded by an obscuring torus. 

We have some significant pieces of information on the 
obscuring materials. By X-ray observations, it is indi- 
cated that most Seyfert 2 nuclei are heavily obscured along 
the line of sight with at least Ay > 10 and sometimes 
A v > 100 (Matt et al. 1996, 1999; Maiolino et al. 1998b; 
Bassani et al. 1999; Risaliti, Maiolino, & Salvati 1999). 
On the other hand, Ay of the nuclear or circumnuclear 
regions is estimated to be between a few and several by 
IR and optical observations (Rix et al. 1990; Roche et al. 
1991; Goodrich, Veilleux, & Hill 1994; McLeod & Rieke 
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1995a; Oliva, Marconi, & Moorwood 1999). Also, it is 
argued that a component of obscuring materials must be 
extended up to > 100 pc in addition to a compact compo- 
nent confined into subparsec scales (Rudy, Cohen, & Ake 
1988; Miller, Goodrich, & Mathews 1991; Goodrich 1995; 
McLeod & Rieke 1995a; Maiolino ct al. 1995; Maiolino & 
Rickc 1995; Malkan, Gorjian, & Tarn 1998). These facts 
may suggest that the distributions of dusty gas around 
an AGN are much more diverse than used to be consid- 
ered, and they have a close relation with circumnuclear 
starburst events. 

Recently, Ohsuga & Umemura (1999, hereafter Paper I) 
have suggested a novel picture for the starburst- AGN con- 
nection, where a large-scale dusty wall of several hundred 
parsecs is built up due to radiation force by a circumnu- 
clear starburst as well as an AGN. This model provides 
a possibility that an AGN type is regulated by circum- 
nuclear starbursts, and also it gives a physical origin of 
extended obscuring materials. However, the dusty wall 
was analyzed in an optically thin regime. Thus, the pos- 
sibility of the strong obscuration with Ay » 1 is still an 
open question in this picture. 

In this paper, by taking the optical depth of dusty gas 
into consideration, we consider the effects of radiation 
force by a circumnuclear starburst as well as an AGN, 
and investigate the stable configuration of dusty gas. To 
begin with, we solve hydrostatic equations coupled with 
the frequency-dependent radiative processes with assum- 
ing spherical symmetry. By this analysis, we demonstrate 
that a geometrically thin, optically thick wall forms due 
to radiation pressure by a starburst. Next, we obtain the 
stable equilibrium configuration of dusty gas walls in two- 
dimensional axisymmetric space. Finally, we attempt to 
classify the AGN type according to the starburst luminos- 
ity and the AGN luminosity. Moreover, taking account of 
the stellar evolution in starburst regions, we discuss the 
time evolution of the circumnuclear structure. In §2, the 
radiation fields and gravitational fields are modeled. In 
§3, the structure of a dusty wall is analyzed in spherically 
symmetric approximation, demonstrating that a geomet- 
rically thin, optically thick wall can form due to radiation 
force by a starburst. In §4, the configuration of dusty 
gas in two-dimensional axisymmetric space is investigated 
with an approximation of geometrically-thin walls, and the 
conditions for the wall formation are given. In §5, based 
upon the present picture, we discuss the implications for 
the starburst-AGN connection. Furthermore, the evolu- 
tion of AGN type which is predicted by the formation of 
the obscuring walls is shown in §6. §7 is devoted to the 
conclusions. 

2. RADIATION FIELDS AND GRAVITATIONAL FIELDS 

Recent observations have been revealed that the circum- 
nuclear starburst regions frequently exhibit ring-like fea- 
tures (Wilson et al. 1991; Forbes et al. 1994; Marconi et 
al. 1994; Mauder et al. 1994; Buta, Purcell, & Crocker 
1995; Barth et al. 1995; Leitherer et al. 1996; Maoz et al. 
1996; Storchi-Bcrgman, Wilson, & Baldmin 1996). There- 
fore, as a component of radiation sources, we assume a 
starburst ring whose bolomctric luminosity, radius, and 
total mass are L$b, Rsb, and M$b, respectively. [Even 
if the starburst regions are localized, the radiation fields 



could be equivalent to that by a ring-like source owing to 
the short rotation timescale (see also Paper I)]. We also 
include an AGN itself, in which the bolometric luminos- 
ity is Lagn and we simply assume an energy spectrum of 
L^ GN oc v- x between 0.01 eV and 100 kcV (Blandford, 
Netzer, & Woltjer 1990). To calculate the bolometric lu- 
minosity and the energy spectrum of the starburst ring, we 
assume the initial mass function (IMF) to be of Salpeter- 
type, 

</>oc K/M )- L35 , (1) 

and the star formation rate (SFR) in the starburst regions 
to be 

SFRcxexpt-^), (2) 

where isB is the elapsed time after the initial starburst. 
The total luminosity is regulated by the mass-luminosity 
relation, (U/Lq) = (m»/M Q ) 3 - 7 , and the mass-age rela- 
tion, t = 1.1 x lO lo yr(m*/M0) -2 - 7 , where m» and Z» are 
respectively stellar mass and stellar luminosity. Here, we 
consider the mass range of 2M© — 40M Q in the IMF, since 
several observations indicate that the IMF is deficient in 
low mass stars with the cutoff of about 2M , and the 
upper mass limit is inferred to be around 40M Q (Doyon, 
Puxley, & Joseph 1992; Chariot et al. 1993; Doane & 
Mathews 1993; Hill et al. 1994; Brandl et al. 1996). 
In addition, we employ the mass-temperature relation, 
(T*/T Q ) = (m»/M Q ) ' 575 , and the energy spectrum of the 
star is assumed to be the blackbody spectrum with the ef- 
fective temperature, T* . We neglect the radiation from the 
supernovae, because the long term-averaged luminosity of 
supernovae is less than 10% of the total stellar luminosity 
of the ring as shown in Paper I. Resultantly, the bolomet- 
ric luminosity of the starburst ring is in proportion to the 
total mass, and decreases with time. 

As for the gravitational fields, we should consider the 
galactic bulge, the central black hole, and the starburst 
ring. But, the gravitational fields in circumnuclear regions 
suffer from some uncertainty. Recent observations provide 
significant information which allows us to model the grav- 
itational fields in circumnuclear regions. For instance, the 
stellar rotation velocity in the Circinus galaxy (a closest 
Seyfert 2 galaxy) shows that the mass distributions in cir- 
cumnuclear regions are not concentrated into a point-like 
object but extended to a few hundred parsecs (Maiolino 
et al. 1998a). Therefore, we consider also such an inner 
bulge-like component. We assume the galactic bulge to be 
an uniform sphere whose mass and radius are Mqb and 
Rqb j the inner bulge also to be an uniform sphere whose 
mass and radius are M\b and -Rib, and the mass of the 
black hole to be Mbh, respectively. 

The observations of IRAS galaxies by Scoville et al. 
(1991) show that the central regions within kilo parsecs 
possess the mass of <1O 1O M . This is comparable to a 
typical mass of galactic bulge. A rotation velocity of a 
starburst ring also implies that the dynamical mass within 
a few hundred parsecs is several 10 9 M Q (Elmouttie et 
al. 1998), which provides an estimation of the mass of 
inner bulge. In addition, Maiolino et al. (1998a) esti- 
mated an upper limit of a putative black hole as several 
10 6 M Q . (We tentatively assume the black hole of 1O 7 M , 
but the overall results are not changed much even if one 
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assumes M B h = 10 6 M Q or 10 s M Q .) By taking these ob- 
servational data into account, we adopt the mass ratio as 
M GB : M m : M BH = 1 : 0.1 : 1CT 3 . On the other hand, 
the mass of the starburst ring, Msb, is supposed to be 
between 1O 7 M and 10 10 M Q . Also, a starforming ring is 
observed around a few hundred parsecs in some galaxies 
(Marconi et al. 1994; Elmouttie et al. 1998). Hence, here 
we assume a starburst ring at a radius of 200 pc. For the 
size ratios, we employ Rgb ■ Rsb ■ Rib =5:1: 0.5. 

3. STRUCTURE OF RADIATIVELY-SUPPORTED OBSCURING 

WALL 

The equilibrium configuration of dusty gas is nearly 
spherically symmetric as shown in the next section. Thus, 
first we investigate the structure of a radiatively-supported 
gas wall inside the starburst ring under an approximation 
of local spherically symmetry. For this purpose, we as- 
sume the gravity by the starburst ring to be constant at 
a part of the wall of interest. Also, the radiation force by 
the starburst ring is assumed to to be a function solely 
of optical depth measured from the outer surface of the 
wall. (Rigidly speaking, both assumptions are verified 
only when the wall is geometrically thin and located at 
the equatorial plane.) Under these assumptions, we solve 
the one-dimensional radiation-hydrostatic equation, cou- 
pled with ionization process and thermal process. 
The radiation-hydrostatic equation is given by 
GM l 1 dP e f X v L$ G "e-T~ SB 
+ J c 47T/2 au + J * 



I 2 



Pg dl 



I grav 



- T »Uv = 0, (3) 



where I is the radius, Mi is the total mass within the ra- 
dius I, p g is the gas density, P g is the gas pressure, \v 
is the mass extinction coefficient of dusty gas, t„ is the 
optical depth of the wall measured from the center, T * otal 
is the total optical depth of the wall, and f g ® v and F^ B 
are respectively the gravity and the radiation flux by the 
starburst ring. It is noted that f g ® v and F^ B provide 
non-spherical components in real situations. The effects 
of the non-sphericity are taken into consideration by eval- 
uating these quantities by two-dimensional axisymmetric 
calculations. 

Recent observations report the existence of a large 
amount of dust, molecular gas, and metals in QSOs (Bar- 
vainis, Antonucci, & Coleman 1992; Ohta et al. 1996; 
Omont et al. 1996; Dietrich & Wilhelm-Erkens 2000). 
Therefore, we suppose the dust-to-gas mass ratio to be 
0.03, which is three times as large as that observed in 
the Solar neighborhood. Since the mass density of dust 
and Thomson scattering are negligible, the extinction is 
given by Xv = (™hXhic™ + ofy/p & . Here, n H is the 
number density of hydrogen nuclei, xhi is the fraction 
of neutral hydrogen, cr™ is the photoionization cross- 
section, and af, is the absorption coefficient of dust. For 
the dust model, we employ the grain size distribution 
as n d (a d ) oc aj 3 ' 5 in a range of [0.01/xm, 1/im] (Mathis, 
Rumpl, & Nordsieck 1977), and absorption cross-section, 
7ra^min[l, (A/27rad)~ 2 ], where ad is the grain radius, and 
the density of solid material within a grain is assumed to 
be l.Ogcm -3 . 

The equation of the ionization balance is 
r 7 Xffl + r ci n H Xffl(l - Xhi) = a re cn H (l - Xhi) 2 , (4) 



where T 7 is the photoionization rate, T" is the collisional 
ionization rate, and a roc is the recombination coefficient. 
T 7 is given by 



p7 



4irl 2 



+ F b 



SB p -(r*' 



1 

hv 



dv, 



with ui, being the Lyman limit frequency, and r cl is 
1.19 x 10- 10 T g 1/2 exp(-1.58 x 10 5 /T g ) (Sherman 1979). 
If a free electron recombines directly to ground state 
of hydrogen, the emitted photon has enough energy to 
cause further photoionization. Since a part of recombi- 
nation photons are absorbed by the dust grain, the value 
a ICC is given by a ICC = a A - (a a - OTsWxHic^/PgXi'L, 
where total recombination coefficient, a a, is well fit- 
ted by a A = 2.1 x 10- 13 (T g /10 4 K)- 1 /2^(i6T g /10 4 K), 
with <p(y) — 0.5(1.7 + lny + 1/6?/) for y > 0.5, or 
y(-0.3 - 1.2 In y) + y 2 (0.5 - lny) for y < 0.5 (Sherman 
1979), and the recombination coefficient to all excited lev- 
els of hydrogen, «b, is approximately related with a a as 
a B = a A exp{-0.487(T g /10 4 K) 1 / 5 }, where T g is the gas 
temperature. 

The equation for the energy balance of the gas is given 

by 

£f + C c ' + Cl e + Cl + Cf, (6) 



nz = ci 



where 7Y 7 is the heating rate of the gas due to direct ra- 
diation by the starburst ring as well as the AGN, and 
£™ c , C c g \ C c g e , C%, and Cf, are respectively the cool- 
ing rates for gas through radiative recombination, free-free 
emission, collisional ionization, collisional excitation, emis- 
sion for metal, and collision with dust. HJ is given by 



K = / XHI< 



f l^ GN C T - SB (r tota._ Ti/ 

\ 47T/2 



(1-^,(7) 



where we neglect the energy transfer by the Compton 
scattering because the timescale is estimated to be ~ 
10 9 yr and it is much longer than the typical lifetime 
of AGN, ~ 10 s yr, in the circumnuclear regions of > 
several 10 pc even in the case of very luminous AGNs 
with ~ 10 13 L Q . £™ c , C c g \ and C c g c are repre- 

sented by £ g oc = 3fcT g a rcc n H (l - Xhi) 2 /2, = 1.42 x 

10- 27 T g 1/2 n H (l - xhi) 2 , = r ci ^ L n H Xffl(l - Xhi), 
£ c g c = 7.5 x 10- 19 exp(-1.18 x 10 5 /T g )n H Xffl(l - X hi) 

(Black 1981), and C\ = 4 x 10- 28 T g 1/2 n H for 10 2 K < 
T g < 10 4 K, or Cl - 10- 22 1 - 120Z T g 30Z n H - 10- 221 n H 
for 10 4 K < T g < 10 5 K (Theis, Burkert, & Hensler 1992), 
where Z is the metal abundance. Here, the metal abun- 
dance is assumed to be 3Z Q , consistently with the sup- 
posed dust-to-gas mass ratio, 0.03. Assuming the colli- 
sional equilibrium of charge on a grain, we can estimate 
Cf = J n d TTa 2 d da d -(7~6 X Bi)(SkT g /TTm p ) 1 / 2 (2kT g -2kT d ) 7 
where m p is the proton mass and Td is the dust temper- 
ature. As for the energy balance for dust, the emission 
from gas is assumed to be absorbed by dust almost on the 
spot. Then, the energy equation for dust is 



HI 



*-d' 



(8) 
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where TC d and £l ec 



arc the heating rates for dust due to the 
direct radiation from a starburst ring as well as an AGN 
and due to the recombination photons from hydrogen, and 
£j is the cooling by dust emission. (This assumption is 
justified by the fact that a resultant gas wall is optically 
thick due to dust opacity.) Til is given by 



1 r ( T AGN.- 



du, (9) 



£ g cc is written as £g CC = (3fcT g /2 + hv^)a lcc n^(l 
Xhi) 2 , and is presented by 



{r(6)C(6)/r(4)C(4)}(fc/^o) 2 T d 6 , where 



= O^A^h) 



r and £ are re- 
spectively the gamma function and Riemann's function 
(Nakamoto & Nakagawa 1994). We note that dust is pri- 
marily heated by the direct radiation, so that the dust tem- 
perature is not affected heavily even if the heating by the 
gas emission is dismissed. Moreover, since the absorption 
cross-section against IR radiation re-emitted from dust is 
much smaller than that against optical or UV radiation, 
we suppose that dust grains work as pure absorbers. 

The above equations and the equation of state for ideal 
gas are solved to give the gas density, the gas temperature, 
the fraction of neutral hydrogen, and the dust tempera- 
ture as functions of radii. The resultant profile of hydro- 
gen number density and Ay measured from the center are 
shown in Figure 1. Here, Lagn = 10 10 L©, M SB = 1O 9 M , 
and Lsb = 8.3 x lO 10 ^© (£sb = 3 x 10 7 yr) are assumed. 
As seen in this figure, there forms an optically-thick wall, 
where the inner surface is pushed outward by the radia- 
tion force by the AGN against the gravity by the inner 
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Fig. 1. — The profile of hydrogen number density and the vi- 
sual extinction Ay measured from the center are shown as func- 
tions of radii, I. Here, Lagn = 1O 1O L , Msb = 10 9 M Q , and 
L S B = 8.3 X 10 10 L Q (t SB = 3 X 10 7 yr) arc assumed. The AGN 
is located at the center, 1 = 0, and the size of the starburst ring 
is assumed to be 200 pc. The interior and exterior surface of the 
wall arc irradiated by the radiation from the AGN and the starburst 
ring, respectively. The density gradient is very steep since the wall 
is compressed due to strong radiation force by the starburst ring as 
well as the AGN. The effective thickness of the wall is less than 10% 
of its radial extension. This figure shows that a geometrically thin, 
optically thick wall forms. 



bulge and the radiation force by the starburst shoves the 
wall inward. The optical depth (Ay) is basically deter- 
mined by the ratio of the net radiation force to the total 
gravity if both AGN and starburst are super-Eddington 
respectively to Mi and to Msb- The density profile of the 
wall has a peak at ~ 50 pc and the density gradients on 
both sides of a peak are very steep. Then, the geometrical 
thickness of the wall is basically determined by the ratio of 
the square of sound speed to the total radiation force. The 
gas temperature is around 10 K due to the photohcating. 
As a result, the effective thickness of the wall is less than 
10% of its radial extension. Hence, it is concluded that a 
radiatively-supported wall can be optically thick and geo- 
metrically thin. If the starburst is more luminous, the wall 
would be thinner because the density gradient would be 
steeper owing to the radiation force. On the other hand, 
if the mass of the starburst ring is less than 10 9 M Q , the 
density profile would be roughly the same, since the force 
fields are mainly determined by the AGN and the inner 
bulge. In this analysis it has turned out that the radia- 
tion force from the AGN as well as the starburst is exerted 
mainly on dust. Therefore, for the mass extinction of the 
dusty gas, we take only the dust opacity into account in 
two-dimensional calculations below. 

Although we have considered the structure of a gas wall 
inside the starburst ring, the structure of an outer wall 
beyond the starburst ring can be understood by a similar 
argument. If an outer wall forms at several hundred par- 
sees, the photoheating does not work effectively to raise 
the temperature ~ 10 4 K, but it becomes around ~ 10 2 K. 
Then, the thickness of the wall is much thinner than the 
inner wall. As for the optical depth (Ay), it is deter- 
mined basically by the ratio of the radiation force to the 
total gravity for super-Eddington luminosity, and there- 
fore could be Ay 3> 1. 

4. CONFIGURATION OF OBSCURING WALLS 

In this section, we examine the equilibrium configura- 
tion of a gas wall and the stability in two-dimensional 
axisymmctric space. We assume the wall to be geomet- 
rically thin as analyzed in the previous section. Under the 
assumption, taking the effects by the optical depth into 
account, we calculate the radiation force and the gravity 
which are exerted on a dusty wall. The radiative flux force 
by the starburst ring is given by 

Xsb Psb 1 - cxp(-t S b/cos6> sb ) 



/sb 



4tt1 2 sb 



-n l dV, (10) 



T S B /cos 6>sb 

at a point of (r, z) in cylindrical coordinates, where i de- 
notes r or z, xsb is the flux mean mass extinction coef- 
ficient averaged over the starburst radiation, tsb is the 
total optical depth of the wall which is estimated with 
using xsb, Psb is the luminosity density in the starburst 
ring, ?sb is the distance from (r, z) to a volume element 
dV of the ring, #sb is the viewing angle from this element, 
and n % is the directional cosine. Similarly, the radiative 
flux force by the AGN is presented by 

fi _ Xagn ^agn 1 - cxp(-t A gn/ cos6> AGN ) i 

J AGN c ^ TAGN / COS $AGN V { ' 

at the same point, where xagn is the flux mean mass ex- 
tinction averaged over the AGN radiation, I is the distance, 
#agn is the viewing angle from the center, and tagn is the 
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total optical depth, which is estimated with using xagn- 
Using equations (10) and (11), the equilibrium between 
the radiation force and the gravity is determined by 



dz 



fsB + /. 



AGN 



(12) 



in the vertical directions (z-directions) and 



dr 



/s: 



SB 



JAGN Jei 



grav 



(13) 



in the radial directions (r-directions), where $ is the ef- 
fective potential, j is the specific angular momentum of 
the dusty gas and /g rav is the gravitational force. Further- 
more, stable configuration should satisfy the conditions 



dz 2 



>0, 



(14) 



and 



dr 2 



> 0. 



(15) 



4.1. Inner Obscuring Wall 

First, we consider the formation of a dusty wall inside 
the starburst ring. The wall is irradiated outward by the 
AGN and inward by the starburst ring. In Figure 2, the re- 
sultant equilibrium branches are shown in the r-z space for 
A v = 3, 5, and 7. Here, Lagn = 10 10 L©, M SB = 10 9 M Q , 
and L SB = 8.3 x 10 10 L o (i SB = 3 x 10 7 yr) are assumed. 
On each curve, the conditions for stable equilibrium in the 
vertical directions [(12) and (14)] are satisfied. However, 
the dashed curves do not satisfy the force balance (13) 
in the radial directions, so that the dusty gas is likely to 
be swung away around the curves because of centrifugal 
force. As a result, it is found that the solid curves are the 
stable equilibrium branches in both vertical and radial di- 
rections, which provide the final configuration of the inner 
obscuring wall. It is noted that the starburst ring plays 
an important role for the force balance and the stability of 
the wall. The AGN as well as the bulge components pro- 
vides only spherically symmetric forces, whereas the radi- 
ation force and the gravity by the starburst ring are not 
spherically symmetric. Therefore, on the tangential plane, 
the azimuthal component of centrifugal force can balance 
only with that of radiation force by the ring (see Figure 
3). To realize this situation, the starburst luminosity is 
required to be super-Eddington, Lsb > 47tcGMsb/xsb- 
Since the effective radiation force by the starburst ring is 
weakened for a large viewing angle 9sb , the radial balance 
breaks down in the vicinity of the rotation axis (z-axis). 



o 
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Fig. 2. — Equilibrium configuration of inner walls is shown in r-z 
space for A v = 3, 5, and 7. Here, L A gn = 10 10 -L Q , M SB = 10 9 Af Q , 
and L S B = 8.3 x 10 10 L o (*SB = 3 x 10 7 yr) are assumed. The solid 
curves represent the stable branches, while the dashed curves are 
radially nonequilibrium branches. The solid curves show the final 
configuration of stable inner walls. Since the effective radiation force 
by the ring is weakened in the vicinity of the rotation axis in the 
case of the inner wall due to a large viewing angle, an opening forms 
and it is shown by a dashed curve. 




AGN 



Starburst 
Ring 

Fig. 3. — Schematic edge-on sketch showing the forces exerted on 
the inner wall as well as the outer wall. Arrows represent the orienta- 
tions of forces, although the lengths of the arrows do not correspond 
to the strength of forces. As a result of numerical calculation, it is 
found that the starburst ring provides nonspherical radiation force 
and gravity, whereas the forces by the AGN as well as the bulge 
components are spherically symmetric. Hence, the walls form if the 
azimuthal (0) component of radiation force by the starburst ring 
balances with that of total force of the gravity by the ring and the 
centrifugal force. 

An opening grows as the optical depth (Ay) increases, be- 
cause of the dilution of the effective radiation force by the 
ring [see equation (10)]. In order that a stable inner wall 
enshrouds the nuclear regions with a large covering factor 
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(opening angle < 90°), Ay is constrained from above as 



XV L SB 

AttcGMsb ' 



(16) 



where xv is the mass extinction at the V band. On the 
other hand, the radiation force by a very luminous AGN 
could blow the dusty gas away from the inner bulge regions 
if the optical depth is not enough to dilute the radiation 
force. This provides a lower bound of Ay as 



Ay > 



XV L AGN 



47tcG(Mib + M BH ) 



(17) 



Because of the time dependence of starburst luminosity, 
the upper bound of Ay given by inequality (16) is a de- 
creasing function of time. For a given model, the allowed 
range of Ay of the inner wall is presented in the Figure 
4 as a function of time, ts&- In the shaded regions, no 



T 



Radially nonequiriblium 




t SB [10 7 yr] 

Fig. 4. — The range of Ay allowed for the inner wall formation. 
In the shaded area, Ay does not satisfy the inequality (16), so that 
the greater part of the wall is out of radial equilibrium. Also, the 
dusty gas is blown away due to strong radiation force by the AGN 
in the hatched area in the case of Lagn = 3 x 1O 1O L0. 

stable equilibrium branches with a large covering factor 
exist and in the hatched regions the radiation force by the 
AGN blows out the dusty gas. If the upper bound in (16) 
is smaller than the lower bound in (17), there is no allowed 
value for Ay . That is to say, no inner wall forms regardless 
of t SB if 

Lagn > (M IB + M BH ) (18) 

Msb 

4.2. Outer Obscuring Wall 

Next, we investigate the formation of an outer wall with 
the radial extension of > i?sB- In Figure 5, the resul- 
tant equilibrium branches of the outer wall are shown in 
the r-z plane. As shown in this figure, a nearly spheri- 
cal wall forms. In this case, the starburst ring as well as 
the AGN irradiate the inside of the wall. Therefore, the 
radiation force by the starburst ring is exerted outward 
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~i 1 1 r 



~i 1 1 r 



L AGN =1.5xlO^L e 



A„=12 




1000 



Fig. 5. — Equilibrium configuration of an outer wall beyond the 
starburst regions is shown in r-z space for two cases of the AGN 
luminosity, say, -Lagn = 3 x 10 11 Lq and 1.5 X 1O 12 L0. Here, 
M SB = 1O 9 M and L S B = 8.3 X IO^Lq are assumed. The solid 
curves represent stable branches, while the dashed curves are radi- 
ally noncquilibrium branches. The outer wall is roughly of spherical 
shape and covers a wide solid angle whenever it forms. 

on the dusty wall, although the inner wall is pushed in- 
ward by the starburst radiation. But, it has been found 
by the present analyses that the radial equilibrium and sta- 
ble condition for the outer wall formation is again given 
by the inequality (16) for a wide range of parameters. It is 
also required that the starburst ring is super-Eddington. 
This can be also understood by the argument of force 
balance in the azimuthal directions at a point on the 
wall. The azimuthal component of the centrifugal force 
can be balanced only by the force by the starburst ring, 
because only the starburst ring provides a spherical force 
(see Figure 3). The upper value of the allowed Ay de- 
creases again with time as the starburst becomes dimmer. 
As a prerequisite condition for the formation of the wall, 
(aagn^agn + Xsb£sb)/47tcG(M S b + M m + M B h) > 1 is 
required, otherwise the dusty gas is not held at the regions 
of > i?iB even in the case of an optically-thin wall. Since 
the radiation force is diluted by the optical depth of the 
wall, the upper bound of Ay is provided as 



Ay < 



Xv (Lagn + £sb) 



(19) 



4ncG(M SB + Mm + M BH ) ' 

for the sustainmcnt of the dusty gas. On the other hand 
the lower bound of Ay is estimated by 

Xv (Lagn + £sb) 



Ay > 



M; 



SB 



Mi 



IB 



BH/ 



(20) 



4ttcG(M gi 

since the radiation force would blow out the dusty gas in 
most regions if the Ay of the wall is smaller than this value. 
Comparing this lower limit to equation (16), we find that 
no outer wall forms if the AGN luminosity satisfies the 
inequality as 

Lagn > (M GB + M m + M BH )^. (21) 

M S B 



Formation of Obscuring Walls by Starburst Radiation 



7 



It is noted that this condition involves the prohibition con- 
dition (18) for an inner wall. Ay is limited from above as 

Xv^sb Xv (^agn + £sb) 



Ay < min 



4ttcGM sb ' 47rcG(M SB + M IB + M BH ) 



(22) 



5. IMPLICATIONS FOR STARBURST- AGN CONNECTION 

From a viewpoint of the formation of radiatively- 
supported obscuring walls, we discuss the connection be- 
tween the AGN type and the starburst events. Based upon 
the conditions of (18) and (21), in Figure 6 the prediction 
for the wall formation is summarized in a diagram of Lagn 
versus Lsb/M$b,9, where Msb,9 is the mass of a starburst 
in units of 1O 9 M . An outer wall forms below a solid line 
and an inner wall as well as an outer wall forms below a 
dashed line. A vertical dot-dashed line shows the bound- 
ary for the formation of the walls. The possible values 
of Ay are also shown by thin dot-dashed lines. Here, it 
is stressed that the possible Ay is practically determined 
by equation (22) for the outer wall. However, if the mass 
supply to an altitude of several 100 pc by a superbubble 
driven by a circumnuclear starburst is insufficient, no outer 
wall of Ay > several might form (see next section). The 
figure shows that for a luminous starburst and a relatively 
faint AGN, both inner and outer walls are built up, and 
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Fig. 6. — The conditions for the formation of obscuring walls 
are shown in a diagram of the AGN luminosity Lagn versus the 
starburst luminosity Lsb, where the starburst luminosity is divided 
by the mass of starburst in units of 10 9 Mq. An outer wall forms 
below a solid line and an inner wall as well as an outer wall forms 
below a dashed line. A vertical dot-dashed line gives the bound- 
ary for the formation of the walls based on the condition that the 
starburst is supcr-Eddington, and the specific values of Ay are also 
shown by thin dot-dashed lines. The figure shows that for a lumi- 
nous starburst and a relatively faint AGN, the nucleus is enshrouded 
by double walls with the total extinction of Ay > 10, and therefore 
it is likely to be identified as a type 2 AGN (Sy 2). If a starburst is 
intrinsically faint or becomes dimmer owing to the stellar evolution, 
neither wall forms and the nucleus is observed as a type 1 AGN (Sy 
1). The age of the starburst based on the present stellar evolution 
model is shown in the upper abscissa. Furthermore, if the AGN is 
much more luminous, the wall formation is prohibited regardless of 
the starburst luminosity (QSO). 



the nucleus could be obscured with the total extinction of 
Ay ~ 10. Then, the nucleus is likely to be identified as a 
type 2 AGN (Sy 2). Also, the outer wall may be consistent 
with the recently observed obscuring material extending 
up to > 100 pc around the nuclei (Rudy, Cohen, & Ake 
1988; Miller, Goodrich, & Mathews 1991; Goodrich 1995; 
McLeod & Rieke 1995a; Maiolino et al. 1995; Maiolino & 
Rieke 1995; Malkan, Gorjian, & Tarn 1998). 

By X-ray observations, the Ay of most Seyfert 2 nu- 
cleus is estimated to be > 10 and sometimes > 100. Such 
a difference between Ay estimated by the X-ray observa- 
tions and Ay assessed by the IR or Optical observations 
might be due to obscuring matter free from dust, located 
just inside the dust sublimation radius of subparsec. The 
structure and formation mechanism of obscuring material 
in subparsec regions is debatable, but it is beyond the 
scope of this paper. 

If a starburst is intrinsically faint, neither an outer wall 
nor inner one is expected to form. Then, the nucleus could 
be observed as a type 1 AGN (Sy 1). These results provide 
a physical explanation for the putative correlation between 
AGN type and host properties whereby Seyfert 2 galaxies 
are more frequently associated with circumnuclear star- 
bursts than type 1 galaxies. In the cases of the moderate 
extinction by the walls may lead to intermediate types of 
AGNs, that is, type 1.9, 1.8, 1.5, or 1.2 according to the 
decrease of the extinction. 

In addition, we have come to a significant conclusion 
that a much higher AGN activity as (21) would preclude 
either wall from forming. Then the nucleus tends to be 
identified as type 1 regardless of the starburst luminosity. 
This result may be closely related to the fact that QSOs 
are mostly observed as type 1 regardless of star-forming 
activity in the host galaxies. 



6. AGN EVOLUTION 

In the previous section, it is shown that the obscuring 
walls form in the case of the luminous starburst and the 
relatively faint AGN. It implies that the obscuration of the 
AGN is time-dependent, since the bolomctric luminosity of 
the starburst ring decreases with time due to stellar evolu- 
tion. In this section, we discuss the evolution of the AGN 
type predicted by this picture. The schematic view of the 
evolution of circumnuclear structure is shown in Figure 7. 

In an early evolutionary stage, a large amount of dusty 
gas could be blown out by superbubbles in the circumnu- 
clear starburst regions. Shapiro & Field (1976), Tomisaka 
& Ikeuchi (1986), and Norman & Ikeuchi (1989) have 
shown that a superbubble of ~ lkpc powered by mas- 
sive OB associations forms. As an observational exam- 
ple, the superbubble of l.lkpc is observed in the immedi- 
ately east of the nucleus of nearby Seyfert 2 /LINER galaxy 
NGC3079 (Veilleux et al. 1994). 

Also, if the mass supply rate given by Norman & Ikeuchi 
(1989) is applied to the case of circumnuclear starbursts of 
Msb ~ 10 9 M Q , the rate into the regions of several 100 pc 
height from the circumnuclear molecular disk due to super- 
bubbles is estimated as several x 10M Q yr _1 . In the case of 
the superbubble of the NGC3079, Duric & Seaquist (1988) 
and Veilleux et al. (1994) evaluated the mass supply rate 
to be >10M Q yr~ 1 and ~ 7M©yr _1 , respectively. As a 
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Fig. 7. — Schematic illustration of the evolution of circumnuclear 
structure. Within ~ 10 7 yr, a large amount of dusty gas is blown 
away due to superbubbles powered by type II supcrnovac in cir- 
cumnuclear starburst regions. The ejected dusty gas is supported 
by the radiation force by the AGN as well as the starburst ring, 
and composes the obscuring walls. Then, the nucleus is obscured 
by both walls and tends to be identified as type 2 AGN (Sy 2). 
Based on the present stellar evolution model, the nucleus tends to 
be observed as type 2 AGN (Sy 2) within several X 10 7 yr, since the 
possible Ay decreases with time due to stellar evolution. Finally, 
the walls disappear and the nucleus are observed as type 1 AGN (Sy 
1). Therefore, the AGN is destined to shift from type 2 to type 1 
in several times 10 7 yr in this picture. However, QSOs are observed 
as type 1 AGN regardless of the age of the starburst, because the 
luminous AGN prevent the formation of the walls. 

result, the dusty gas of ^>1O 8 M would be supplied within 

a typical lifetime of a superbubble (~ 10 7 yr) and com- 
poses the obscuring walls in the next stage. However, an 
outer obscuring wall with Ay > 10 may not form, because 

it requires the dusty gas of >1O 8 M . On the other hand, 

the required mass for an inner wall is merely ~ 10 6 Mq. 
The formation timescales of an inner wall and that of an 
outer wall are several x 10 5 yr and several x 10 6 yr, respec- 
tively. Both timescales are shorter than the evolutionary 
timescale of the starburst on the order of ~ 10 7 yr. Thus, 
the nucleus is likely to be obscured by the inner wall with 
Ay > several as well as the outer wall with Ay < a few, 
whereas the possible Ay decreases with time due to the 
stellar evolution in the starburst regions. Based on the 
present stellar evolution model, the nucleus is heavily ob- 
scured and tends to be observed as type 2 AGN (Sy 2) 
within several xlO 7 yr. Also, the inner wall would not 
be destroyed even if the superbubbles occur in this stage, 
since the bubble diameter on the equatorial plane is com- 
parable to the scale height of the molecular disk (Schiano 
1985; Mac Low & McCray 1988; Tomisaka & Ikeuchi 1988; 
Mac Low, McCray, & Norman 1989). 

In the final stage, the obscuring walls disappear due to 
the diminution of the bolometric luminosity of the circum- 
nuclear starburst. Then, the nucleus would be identified 
as type 1 AGN (Sy 1). To conclude, the AGN type is time- 
dependent and is destined to shift from type 2 to type 1 in 
several times 10 7 yr. It is stressed that no obscuring wall 
forms in the case of the luminous AGN like QSOs so that 
type of QSO does not vary, and which is observed as type 
1 regardless of the age of the starburst. 

Here we have assumed that AGN activity and the cir- 



Umemura 

cumnuclear starbursts are simultaneous events. From an 
observational point of view, it seems to be reasonable, be- 
cause it is revealed that AGN events are frequently ac- 
companied by starbursts in a variety of data (Scoville et 
al. 1986; Soifer et al. 1986; Heckman et al. 1989, 1995; 
Marconi et al. 1994; Neff et al. 1994; Cid Fernandes 
& Terlevich 1995; Genzel et al. 1995; Maiolino et al. 
1995, 1998a; Oliva et al. 1995; Storchi-Bergmann et al. 
1996; Rodriguez-Espinosa, Rudy & Jones 1987; Elmouttie 
et al. 1998; Cid Fernandes, Storchi-Bergmann, & Schmitt 
1998). In addition, as a physical solution which links the 
two events, Norman & Scoville (1988) have suggested a 
model, whereby the interstellar medium which is released 
by post-main-sequence stars in starburst regions feeds a 
central black hole. Also, the radiatively driven mass ac- 
cretion onto a black hole due to the radiation drag, which is 
proposed by Umemura, Fukue, & Mineshige (1997, 1998), 
Fukue, Umemura, & Mineshige (1997), and Ohsuga et al. 
(1999), would link the two events. 

In the present analysis, we have assumed the interstel- 
lar gas containing a large amount of dust which is inferred 
in QSOs (Barvainis, Antonucci, & Coleman 1992; Ohta et 
al. 1996; Omont et al. 1996; Dietrich & Wilhelm-Erkens 
2000). However, the dust grains could be destroyed in 
shock-heated hot gas and they could be reproduced in 
cooled superbubble winds. The dust properties around 
AGNs are under debate also from an observational point 
of view. The destruction, formation, and growth of the 
dust grains would make obviously significant effects not 
only in the present picture but also in a generic problem 
of obscuration. Thus, this issue should be considered more 
carefully in the future analysis. 

Finally, we discuss the stabilities of the obscuring walls. 
First of all, the walls are unlikely to be subject to the 
Ray leigh- Taylor instability. As shown in Figure 1, the den- 
sity gradient of the inner wall is positive inside the density 
peak and negative outside the peak. As shown in Figure 3, 
the acceleration works in the same directions as the density 
gradient in both sides. Thus, the Rayleigh- Taylor instabil- 
ity does not occur at the inner wall. In the case of the outer 
wall, the radiation force by the AGN and the starburst 
ring, which pushes up the wall, is the dominant force at the 
inner surface, but the dusty gas is attracted by the gravity 
at the outer surface. As a result, both walls turn out to be 
stable to the Rayleigh- Taylor mode. However, the walls 
might be subject to the other instabilities. When the den- 
sity perturbations occur in the walls, the thick parts of the 
walls go down by the gravity, whereas the radiation force 
lifts the less thick parts, since the radiation force is ineffec- 
tive for heavily optically-thick parts. Then, the walls are 
likely to have more complicated configurations. Moreover, 
thermal instabilities might break out, because the radia- 
tive cooling works at the dense regions effectively, whereas 
the gas remains as warm as ~ 10 4 K in the optically thin 
regions due to radiative heating. Simultaneously, the self- 
gravitating instability might grow at the cool and dense 
parts of the walls. Also, the actual circumnuclear star- 
burst is not perfectly ring-shaped but more or less clumpy 
and therefore makes non-axisymmetric radiation fields, al- 
though the wheel effects by the rotation could smear out 
non-axisymmetric effects to some degree (see also Ohsuga 
& Umemura 1999). In addition, the dusty walls are pos- 
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sibly damaged due to the collision of starburst winds or 
dusty gas expelled from mass losing stars inside the walls. 
Consequently, the walls could be distorted and fragment 
into small pieces. How the gas clouds resulting from the 
fragmentation obscure the AGN is an interesting issue, but 
the details are not clear before radiation-hydrodynamical 
simulations are performed. 

7. CONCLUSIONS 

In the present analysis, by taking optical depth due to 
dust opacity into consideration, we have analyzed the sta- 
ble equilibrium configuration of dusty gas in the circum- 
nuclear regions in radiation fields by a starburst and an 
AGN. It has been found that the radiation pressure by a 
circumnuclear starburst sustains an inner obscuring wall 
of several ten parsecs and an outer obscuring wall of sev- 
eral hundred parsecs. The total extinction of the obscuring 
walls has turned out to be around ten. Then, the nucleus 
would be identified as type 2. If a starburst is intrinsically 
faint or becomes dimmer owing to the stellar evolution 
there, neither an outer wall nor inner one would not form. 
Then, the nucleus could be observed as a type 1 AGN (Sy 
1). The results on outer and inner walls share a certain 
similarity in that the strong radiation force by the AGN 



prevents the formation of an outer wall as well as an in- 
ner wall. These results provide a physical explanation for 
the correlation between AGN type and host properties, 
whereby Seyfert type 2 galaxies are more frequently asso- 
ciated with circumnuclear starbursts than type 1 galaxies, 
whereas QSOs arc mostly observed as type 1 even though 
vigorous star-forming activities are often observed in the 
host galaxies. Also, the large-scale outer wall of several 
hundred parsecs is consistent with the observations on ob- 
scuring material extending up to > 100 pc around the 
nuclei. Finally, this model predicts the time evolution of 
the AGN type from type 2 to type 1 in several times 10 7 
yr as a circumnuclear starburst becomes dimmer due to 
stellar evolution. 
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